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Abstract —  The  transformer  temperature  is  one  of  the 
main  variables  of  interest  in  its  manufacture  and 
operation,  since  it  interferes  in  the  useful  life  of  the 
transformer.  In  this  sense,  the  present  work  proposes  a 
thermal  simulation  for  the  temperature  estimate  of  a  dry- 
type  transformer.  Initially  a  thermal  analysis  was 
performed  from  experimental  measurements  of 
temperatures  of  a  dry -type  power  transformer  of500kVA 
for  different  conditions  of  load.  Posteriorly,  a  thermal 
simulation  was  proposed  using  finite  element  theory. 
Thus,  the  heat  diffusion  equation  was  used,  with  the 
following  boundary  conditions:  convection  and  radiation 
equations;  characteristics  of  the  materials  used;  the 
measurement  data  and  the  dimensions  of  the  transformer. 
FEMM  2D  software  was  used  for  the  proposed 
simulation.  Finally,  in  order  to  validate  the  proposed 
analyzes,  experimental  measurements  were  compared 
with  the  values  obtained  in  the  thermal  simulation.  The 
results  of  the  thermal  simulation  showed  agreement  with 
the  experimentally  measured  values. 

Keywords —  Dry-type  transformers,  Experimental 
measurement,  Finite  elements,  Thermal  simulation. 

I.  INTRODUCTION 

Dry -type  transformers  are  widely  used  in  electric  power 
distribution  systems  and  are  recommended  for  smaller 
indoor  installations  that  require  safety  and  reliability, 
require  less  maintenance  and  less  damage  to  the 
environment  [1]. 

The  design  and  overload  capacity  of  transformers  are 
strongly  influenced  by  thermal  performance  and  behavior. 
The  thermal  influence  on  the  power  transformers  is  very 
relevant,  and  one  of  the  main  interest  factors  in  the 
operation  of  transformers  is  the  temperature,  since  it 
interferes  with  the  aging  of  the  insulation  of  the  winding 
and,  consequently,  its  useful  life.  Therefore,  the 
temperature  monitoring  is  essential  to  evaluate,  in 
general,  the  wear  of  the  transformer  [2]. 

Studies  have  shown  that  the  temperature  is  one  of  the 
major  causes  of  deterioration  of  insulation  material  of 
power  transformers  and  thermal  stress.  This  thermal 
stress,  caused  by  heat,  is  a  condition  where  the  insulation 


material  is  affected  by  the  temperature  of  the  environment 
in  which  it  is  located,  when  this  temperature  is  at  extreme 
levels,  and  which  has  resulted  in  electrical  faults  in  the 
distribution  systems.  A  temperature  distribution  and 
overload  or  overheating  is  a  particular  interest  of 
manufacturers  and  customers  with  the  purpose  of 
prolonging  the  life  expectancy  of  transformers  [3],  [4]. 

Also,  in  the  design  stage,  in  order  to  estimate  the 
temperature,  there  are  several  numerical  methods  of 
thermal  simulation,  capable  of  dealing  with  complex 
geometries  and  that  become  important  tools  for  the 
transformer  designer.  In  this  sense  the  finite  element 
method  (FEM)  comes  being  widely  applied  after  the 
evolution  of  the  computational  systems  [5]. 

Considering  FEM  this  work  presents  an  analysis 
through  the  simulation  of  the  results  of  theoretical  and 
experimental  investigations  of  the  thermal  behavior  of  a 
dry -type  power  transformer,  when  this  is  subject  to 
different  load  conditions. 

These  investigations,  both  theoretical  and 
experimental,  take  into  account  several  elements  such  as: 
transformer  specification,  thermal  characteristics  of  the 
materials  involved,  transformer  construction  aspects;  as 
well  as  the  data  collected  in  the  experimental 
measurement,  including  the  temperature  at  several  points 
of  the  transformer,  applied  power  and  ambient 
temperature. 

So,  it  is  possible  to  perform  of  dry -type  power 
transformer  simulations  to  estimate  its  temperatures  at 
various  points  and  obtain  thermal  images  of  the  variation 
and  distribution  of  temperature  that  occurs  internally  and 
externally  to  the  equipment.  In  this  way,  the  data 
collected  in  the  experimental  measurement  can  be 
compared  with  the  results  obtained  in  simulations. 

II.  BACKGROUND 

Some  analytical  and  experimental  studies  of  the 
temperature  distribution  in  some  types  of  dry -type 
transformers  have  been  presented  in  the  literature. 

To  this  end  a  thermal  mathematical  model  that 
calculates  of  the  hottest  spot  temperature  and  its  location 
in  dry-type  transformer  windings  was  performed  by  [6]. 
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This  approach  was  based  on  three  mechanisms  of  heat 
transfer:  conduction,  convection  and  radiation.  However, 
the  finite  difference  method  was  applied,  and  circular  or 
oval  forms  of  transformer  geometry  suggested  finite 
element  analysis  to  improve  numerical  accuracy.  The  heat 
conduction  equation  was  also  solved  by  the  finite 
difference  method  by  [7]  in  order  to  obtain  the 
temperature  distribution  for  the  dry -type  transformer  at 
steady  state. 

The  two-dimensional  (2D)  temperature  distribution 
was  calculated  in  [8]  using  finite  element  analysis, 
addressed  by  CFD  (Computational  Fluid  Dynamics) 
analysis  and  using  commercial  ANSYS  software.  This 
work  did  not  consider  heat  dissipation  by  thermal 
radiation. 

Thermal  analysis  and  temperature  distribution  of  the 
transformers  were  also  performed  by  [9]  using  a 
mathematical  model  using  the  finite  difference  method 
and  cylindrical  coordinates.  To  validate  the  model,  the 
results  were  compared  with  the  experimental  data 
measured  from  an  800kVA  transformer. 

Another  2D  thermal  study  of  the  cast-resin  dry-type 
transformer  was  solved  in  [10]  using  cylindrical 
coordinates  and  the  FEM  to  obtain  temperatures  in  any 
location  of  the  transformer.  Basic  modes  of  heat  transfer, 
such  as  conduction,  convection,  and  radiation  were  also 
used  to  obtain  the  steady-state  temperature  distribution  of 
the  transformer. 

More  recently,  a  coupled  thermal-electromagnetic 
model  for  disc-type  and  foil-type  winding  to  propose 
temperature  distributions  using  FEM  was  performed  by 
[11]. 

The  proposed  model  was  based  on  mechanisms  of 
heat  transfer,  i.e.  conduction  heat,  convection  and 
radiation,  in  addition  to  some  electromagnetic 
mechanisms.  All  the  cooling  surfaces  were  identified  and 
heat  transfer  coefficients  for  each  surface  were  presented. 
The  model  was  applied  to  a  2,000kVA  dry -type  power 
transformer  under  different  load  conditions.  Experimental 
temperatures  were  measured  with  thermocouples  and 
infrared  thermometers  and  used  to  compare  the  results  of 
finite  element  simulations.  In  addition,  in  this  model, 
commercial  COMSOL  software  was  used  for  2D  thermal 
modeling  [11].  The  results  of  [6] -[11]  showed  reasonable 
agreement  between  the  calculated  data  and  the 
experimental  data. 

In  this  sense,  the  present  paper  presents  a  thermal 
simulation  of  500kVA  dry -type  transformer  at  steady 
state,  with  an  approach  based  on  the  three  mechanisms  of 
heat  transfer:  conduction,  convection  and  radiation;  in 
cartesian  coordinates,  using  the  MEF  as  numerical 
method  and  the  FEMM  (Finite  Element  Method 


Magnetics)  2D  software  to  perform  the  thermal 
simulation. 


III.  THERMAL  STUDY  IN  DRY-TYPE 
TRANSFORMERS 

The  general  equation  of  heat  diffusion  in  cartesian 
coordinates  2D  and  steady  state  is  represented  in  (1),  in 
which  it  provides  a  basic  tool  for  the  analysis  of  heat 
conduction.  From  the  solution  of  this  equation  it  is 
possible  obtain  the  temperature  distribution  T  (x,  y)  of  a 
steady-state  material  [12]. 


d  /  dT\  d  /  dT\ 

ai(k«a^  +  ap(lt-a^  +  't  =  0 


(1) 


Where  kx  and  ky  are  the  thermal  conductivity 
[W/(mK)]  in  the  directions  of  x  and  y,  q  is  the  heat 
generation  by  volume  [W/m3]  and  T  is  the  temperature 
[K]. 

According  to  [13],  the  term  heat  generation  ( q )  is 

defined  in  (2)  as  the  amount  of  thermal  energy  that  is 
generated  per  unit  volume,  which  in  turn  is  determined  by 
the  amount  of  power  loss  in  a  unit  volume. 


R  = 


(2) 


Where  Eger  is  the  total  rate  of  heat  generation  [W],  1 

is  the  electric  current  [A],  R  is  the  electric  resistance  [Q], 
and  Vis  the  conductor  volume  [m3]. 

The  total  rate  of  heat  generation  is  a  function  of  the 
amount  of  current  flow  in  a  conductor  and  the  resistance 
of  said  conductor. 

The  conductor  resistance  is  calculated  according  to 
(3). 


Where  pe  is  the  electrical  resistivity  [f2.m],  /  is  the 
mean  conductor  length  [m]  and  Ac  is  the  conductor  area 
[m2]. 


3.1  Boundary  Conditions 

When  no  boundary  conditions  are  explicitly  defined,  each 
boundary  imposes  an  isolated  condition.  However,  a  non- 
derived  boundary  condition  must  be  defined  somewhere 
(or  the  potential  must  be  defined  at  a  reference  point  in 
the  domain)  so  that  the  problem  has  a  unique  solution 
[14].  The  convection  (4)  and  radiation  (5)  equations  were 
used  as  boundary  conditions  of  (1). 
dT 

k~+h(T-T0)  =0  (4) 
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k  —  +  £CJ(T4-  T04)  =  0  (5) 

on 

Where  h  is  the  convective  heat  transfer  coefficient 
[W/(rtf.K)],  To  is  the  ambient  temperature  [K],  s  is  the 
emissivity  of  the  surface,  o  is  the  Stefan -Boltzmann 
constant  [W/(n^.K4)]  and  n  represents  the  direction 
normal  to  the  boundary. 

IV.  EXPERIMENTAL  MEASURES 
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The  transformer  used  for  the  experimental  measurements 
was  500kVA  from  the  company  Comtrafo,  which  has  a 
transformation  ratio  of  13.8kV  to  220/127V,  and  it  is  in 
full  use  on  the  Sao  Bernardo  do  Campo  campus  of  the 
Federal  University  of  ABC,  as  ilus traded  in  Fig.  1. 


Fig.  1.  Dry-type  transformer  used  for  the  experimental 
measurements.  Source :  Authors. 

The  equipment  used  for  temperature  measurements 
was  a  thermal  imager,  also  known  as  a  thermal  imaging 
camera  and  a  temperature  sensor  of  type  PT100  installed 
in  the  low  voltage  winding  of  the  transformer. 

A  three-phase  power  analyzer  was  used  to  measure 
electrical  quantities,  such  as:  voltage,  current  and  power. 
The  quantities  were  collected  at  the  same  time  that  the 
temperature  data  was  measured. 

Fig.  2  illustrates  some  of  the  equipment  used  in  the 
experimental  measurements. 


Fig.  2.  Equipment  used  for  the  reading  of  temperature 
(a)  and  electrical  quantities  (b).  Source:  Fluke  [15]. 

For  the  interpretation  and  treatment  of  the  collected 
data,  the  specifical  software  analyzer  was  used. 

Since  the  thermal  behavior  is  analogous  in  the  three 
phases  (coils)  of  the  transformer,  both  the  measurements 
and  the  simulations  were  performed  in  only  one  of  the 
phases. 

Fig.  3a  [16]  illustrates  the  details  of  one  of  the 
transformer  phases  used  for  collecting  measurements  and 
Fig.  3b  shows  the  thermal  variation  found  on  the  external 
surface  using  the  thermal  imager. 


Fig.  3.  Detail  of  one  of  the  phases  of  the  transformer 
used  for  measurement  (a)  and  its  thermal  variation  for  a 
given  power  measured  in  °C  (b).  Source:  (a)  Comtrafo 
[16],  (b)  Authors. 

Since  the  thermal  imager  has  a  limited  visible  screen 
and  there  were  space  constraints  at  the  location  where  the 
transformer  was  located,  the  union  of  four  image 
segments  of  the  thermal  imager  was  used,  all  being  of  the 
same  scale,  thus  forming  a  single  image  of  thermal 
variation  on  the  surface  transformer,  as  shown  in  Fig.  3b. 

TABLE  1  shows  the  data  collected  from  the  electrical 
energy  analyzer,  for  four  different  charging  situations 
(powers). 
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Table.  1:  Data  collected  from  the  energy  analyzer  for 


different  loading  potences. 

Reading 

Voltage 

Current 

Power 

[V] 

[A] 

[W] 

1 

128.8 

63.66 

8,200 

2 

128.4 

74.83 

9,608 

3 

128.7 

97.44 

12,540 

4 

128.5 

130.82 

16,810 

Source:  Authors. 

TABLE  2  shows  the  temperature  results  collected 
from  the  thermal  imager  and  the  thermal  sensor 
considering  the  four  different  loading  powers  of  TABLE 
1. 

Table.2:  Data  collected  from  temperature  meters  for 


different  loading  potences. 


Reading 

Ambient 

Temp. 

[°C] 

Thermal 

Imager  [°C] 

Thermal 

Sensor 

[°C] 

Lower 

Upper 

1 

23 

31.3 

33.8 

50.3 

2 

27 

31.7 

34.5 

50.6 

3 

31 

32.0 

34.8 

50.7 

4 

28 

31.8 

34.9 

50.8 

Source:  Authors. 


V.  COMPUTATIONAL  SIMULATION 

As  a  pre-processing  stage  and  for  the  realization  of  the  2D 
design  of  the  transformer  in  the  FEMM  software,  the 
constructive  characteristics  provided  by  the  manufacturer 
and  the  constant  dimensions  proposed  in  [9]  were  used. 

The  creation  of  core  geometry  and  low  and  high 
voltage  windings  (LV  and  HV)  were  performed  only  in 
one  of  the  phases  of  the  transformer. 

Considering  the  symmetry  of  the  sides  of  the  coil  and 
to  reduce  the  speed  of  the  system  solution  process,  the 
simulation  was  performed  in  only  one  side  of  the  coil,  as 
ilustrated  in  Fig.  4. 

The  data  of  the  thermal  characteristics  of  the  materials 
used  in  the  core,  low/high  voltage  windings  and  in  the 
insulation  material,  such  as:  M4  silicon  steel,  aluminum 
alloy  1,350  and  epoxy  resin,  respectively,  and  their 
variations  of  thermal  conductivity  in  relation  to  the 
increase  in  temperature  were  also  provided  by  the 
manufacturer.  These  thermal  characteristics  of  the 
materials  used  and  the  ambient  air  were  considered  in  the 
simulations. 

As  boundary  conditions,  the  variables  used  in  solving 
the  Equation  (4)  were:  convective  coefficients  and 
ambient  temperature.  For  the  resolution  of  (5)  the 
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variables  involved  were:  emissivity  and  ambient 
temperature. 

cor# 


1  LV  wi  n  d  i  no 


Fig.  4.  2D  drawing  of  one  side  of  the  transformer  coil  in 
the  FEMM  software.  Source:  Authors. 

In  order  to  meet  the  required  variables  in  (4),  the 
surface  where  convective  heat  transfer  occurs  was 
identified  and  divided  into  sections  along  the  height,  these 
partitions  being  uniform  throughout  the  low  and  high 
voltage  winding.  The  values  of  the  upper  and  bottom 
convective  coefficients  of  [11]  were  used  as  reference  for 
the  simulation.  The  convective  coefficients  of  the  interval 
between  the  upper  and  bottom  limits  were  calculated 
based  on  the  temperature  difference  and  the  number  of 
sections  considered. 

In  (5),  for  the  emissivity  values,  the  data  of  [13]  was 
used  as  reference  and  approximation  of  the  values.  The 
emissivity  values  are  based  on  the  temperature  of  300K. 
The  values  of  the  ambient  temperatures  for  both  (4)  and 
(5)  were  presented  in  TABLE  2. 

The  calculated  values  of  the  total  rate  of  heat 
generation  [W]  were  also  considered  in  the  FEMM 
software.  These  values  are  calculated  according  to  (2),  in 
which  they  vary  with  current  and  electrical  resistance. 
The  electric  current  (7)  was  collected  directly  in  the 
transformer  through  the  electric  energy  analyzer,  as 
verified  in  TABLE  1.  The  equation  (3)  was  used  to 
calculate  the  resistance  of  the  conductor. 

To  finalize  the  pre-processing,  the  computational 
domain  was  adapted  and  the  mesh  refinement  was  carried 
out,  which  boils  down  to  increasing  the  minimum  alio  wed 
angle  for  small  triangles  to  be  created.  Thus,  the  meshes 
were  created,  and  the  domains  of  the  problem  were 
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divided  into  finite  elements  using  Lagrange  triangular 
elements,  which  is  the  FEMM  standard  element  type. 

Fig.  5  shows  the  meshes  of  all  domains,  including 
core  and  low  and  high  voltage  windings. 

The  processing  step  consists  in  solving  the  problem, 
based  on  the  parameters  and  variables  defined  in  the  pre¬ 
processing.  The  software  performs  internally  the  matrix 
representation  of  each  element,  formation  of  the  global 
coefficient  matrix,  application  of  the  boundary  conditions 
and,  finally,  the  solution  of  the  system 
After  this,  at  the  final  step,  post-processing,  the  graphical 
and  numerical  results  are  presented. 


Fig.  5.  Finite  element  domains  divided  into  triangular 
elements.  Source:  Authors. 


VL  RESULTS  AND  DISCUSSIONS 

The  thermal  simulation  results  obtained  in  the  post¬ 
processing  step  are  visualized  in  Fig.  6. 

In  the  highlighted  approach  of  Fig.  6  it  is  seen  that  the 
temperatures  in  the  upper  part  are  higher  than  in  the 
bottom  part,  this  is  because  different  convective 
coefficients  were  applied  in  the  sections  of  the  windings 
surfaces,  being  in  agreement  with  the  variations  of  the 
values  collected  in  the  experimental  measurement  for 
high  voltage  windings. 

The  temperature  data  obtained  in  Fig .  6  are  in  Kelvin  and 
the  TABLE  3  presents  the  values  of  the  points  in  °C. 


3.329**002  :  >3.346**002 
3.312**002  :  3.329e*002 
3.294e*002  :  3.312**002 
3.277**002  :  3.294«*002 
3.260**002  :  3.277**002 
3.242e*002  :  3.260**002 
3. 225**002  :  3.242*+002 


3.208**002  : 
3. 190* *002  : 
3.173**002  : 
3.156**002  : 
3.139**002  : 
3.121**002  : 
3.104**002  : 
3.087e*002  : 
3.069**002  : 


3.225**002 
3.208**002 
3.190**002 
3.173**002 
3.156**002 
3.139**002 
3.121**002 
3.104*  *002 
3.087**002 


3.052**002  :  3.069**002 
3.035**002  :  3.0S2e*002 
3.017e*002  :  3.035e*002 
<3.000e*002  :  3.017**002 


Density  Plot:  Temperature  (K) 


Fig.  6.  Visualization  of  thermal  simulation  results. 
Source:  Authors. 


The  TABLE  3  also  presents  the  comparison  between 
the  temperatures  measured  experimentally  with  the  values 
obtained  in  the  thermal  simulation,  as  well  as  the  absolute 
and  relative  error.  From  TABLE  3  it  is  possible  to 
observe  low  absolute  and  relative  errors  and  that  there 
was  a  direct  relationship  of  the  transformer  operating 
temperature  increase  when  compared  to  the  transformer 
power  increase  and  the  ambient  temperature. 


Table  3.  Comparative  of  temperature  measures  and 
simulated  and  their  respective  errors. 


Pow 

Temperature 

er 

[W] 

Site 

Measure 

d[°C] 

Simulat 

ed  [°C] 

Absolut 

e  Error 

[°C] 

Relativ 

e  Error 

Surface 

HV 

31.30 

33.75 

2.45 

8% 

Bottom 

8,2 

Surface 

HV 

33.80 

36.85 

3.05 

9% 

Upper 

LV 

Winding 

50.30 

54.90 

4.60 

9% 

Surface 

HV 

31.70 

34.70 

3.00 

9% 

9,60 

Bottom 

8 

Surface 

HV 

34.50 

37.65 

3.15 

9% 

Upper 

www.iiaers.com 


Page  |  163 


International  Journal  of  Advanced  Engineering  Research  and  Science  (IJAERS)  [Vol-6,  Issue-3,  Mar-  2019] 

https://dx.doi.or a/ 1 0.221 61/iiaers.6.3.1 9  ISSN:  2 349 -649 5 (?)  \  2456-1908(0) 


LV 

Winding 

50.60 

54.70 

4.10 

8% 

Surface 

HV 

32.00 

34.80 

2.80 

9% 

Bottom 

12,5 

Surface 

4 

HV 

34.80 

37.85 

3.05 

9% 

Upper 

LV 

Winding 

50.70 

55.50 

4.80 

9% 

Surface 

HV 

31.80 

34.50 

2.70 

8% 

Bottom 

16,8 

Surface 

1 

HV 

34.90 

37.85 

2.95 

8% 

Upper 

LV 

Winding 

50.80 

55.60 

4.80 

9% 

Source:  Authors. 


The  Fig.  7  illustrates  the  data  presented  in  TABLE  3, 
which  shows  the  measured  and  simulated  temperature 
variations  for  the  power  of  8,200W. 

♦  Measured  ■  Simulated 
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Fig.  7.  Measured  and  simulated  temperatures  with  load 
power  of 8,200W.  Source:  Authors. 

Comparing  Fig.  3b,  image  obtained  experimentally, 
with  Fig.  6,  simulated  image,  it  is  verified  that  both  have 
a  vertical  temperature  gradient.  This  characteristic  occurs 
because  the  convective  heat  transfer  is  larger  in  the 
bottom  part  and  smaller  in  the  upper  part,  obtaining 
higher  temperatures  on  the  upper  surface  and  lower  on  the 
bottom  surface.  Further,  it  becomes  apparent  that  the  air 
temperature  is  heated  with  increasing  axial  length, 
whereby  the  natural  convection  cooling  effect  is  therefore 
greater  in  the  bottom  windings. 

It  can  also  be  seen  in  Fig.  7  that  the  values  measured 
experimentally  compared  to  the  simulated  values  for  the 


same  power  also  showed  agreement  on  the  temperature 
variation  on  the  surface  and  on  the  transformer  windings. 

6.1  Error  Analisis 

The  errors  found  and  described  in  TABLE  3  may  be 
related  to  the  conditions  considered  in  the  thermal 
simulation  or  even  by  inaccuracies  in  the  experimental 
measurements. 

Thus,  effects  of  some  factors  on  the  temperature 
distribution  of  the  measurements  collected  experimentally 
and  the  simulations  carried  out  in  FEMM  software  may 
have  influenced  the  results,  such  as:  inaccuracies  and 
uncertainties  involving  the  temperature  readings  by  the 
thermal  imager  and  the  sensor,  in  the  experimental  part; 
and  errors  caused  by  numerical  calculations  due  to  some 
computational  operations,  such  as  rounding  and 
estimation,  in  the  theoretical  part.  Other  prominent  factors 
were: 

I.  Some  measurements  of  the  size  /  dimensions  for 
the  formation  of  the  transformer  geometry  were 
approximated,  also  contributing  to  the  divergence  of 
measured  and  simulated  values; 

II.  The  simulation  was  performed  in  the  2D  domain 
and  in  cartesian  coordinates,  and  the  3D  computational 
domain  and  cylindrical  coordinates  could  present  more 
accurate  values,  since  they  are  closer  to  the  actual 
configurations  and  geometries  of  the  transformer; 

IE.  Only  the  thermal  study  of  the  transformer  was 
considered,  not  taking  into  account  the  electromagnetic 
interferences  in  the  winding  and  the  core  of  the 
transformer; 

IV.  The  values  of  the  boundary  conditions  for  the 
convective  coefficients  and  emissivities  were  also 
approximated  for  the  geometry  and  applied  materials, 
respectively; 

V.  The  conductor  resistivity  variation  was  not  taken 
into  account  with  increasing  winding  temperature  and 
ambient  temperature,  considering  that  the  resistivity,  and 
consequently  the  resistance,  vary  with  the  increase  in 
temperature  over  time. 

VIL  CONCLUSION 

In  this  work  it  was  possible  to  demonstrate  that  the  finite 
element  analysis  provides  an  efficient  tool  for 
determining  the  temperatures  along  a  dry-type 
transformer  under  load.  Furthermore,  the  heat  transfer 
analysis  allows  to  predict  the  temperature  in  several 
locations  of  the  transformer  and  the  temperature 
distribution  considering  different  operating  conditions. 

Thus,  considering  that  the  transformer  was  simulated 
in  two-dimensional  geometry,  with  some  characteristics 
and  measures  of  approximate  size/dimensions,  it  has  been 
that  the  thermal  simulation  showed  agreement  with  the 
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temperature  variations  of  the  transformer  measured 
experimentally. 

Finally,  the  results  of  this  study  may  provide 
manufacturers  with  a  better  understanding  of  the  nature  of 
heat  transfer  in  transformers,  and  may  help  to  deal  with 
thermal  stresses  more  effectively. 
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